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Abstract. - The multipole expansion technique is applied to one of the largest magnetic 
molecules, Peso. The molecule's dipole, toroid and quadrupole magnetic moments are equal 
to zero (in the absence of magnetic field) so the multipole expansion starts from the octopole 
moment. Probably the Peso molecule is the most symmetrical magnetic body synthesized 
so far. The magnetization process is considered theoretically in different geometries. Some 
components of the octopole moment experience a jump while the magnetization rises linearly 
up to its saturation value. An elementary octopole moment consisting of four magnetic dipoles 
is proposed as a hint for designing of an experiment for measurement of octopole magnetic 
moment components. 



Introduction. - Multipole expansion of the electrical field generated by a number of 
electrical charges has been widely used in theoretical physics. However, multipole expansion 
of magnetic fields generated by a system of permanent currents is less known. The starting 
point of such an expansion is a well known formula for the vector potential of the magnetic 
field: 

i.iv^i^, (1) 

^r\R-ri\ 

where Ci and vl stand for a charge and velocity of the i-th charge. The vector potential A 
can be expanded in a power series of the |7^/|i?| ratio if an interesting point is rather distant 
from the system of charges, i.e. \r\/\R\ <^ 1. The first term of the expansion in question is 
A' = [MR]/R'^ where M = J2i ^d'^if^] = 5J J2i is the magnetic dipole moment of the 
system, U are orbital moments of the electrons/ions and are their masses. The spin si 
degrees of freedom can be taken into account by the substitution h U + gisl. In this paper 
we will be concerned with the spin magnetism only. A distinctive feature of the magnetic field 
potential's ^ multipole expansion is that it is subdivided into two families: magnetic and 
toroid moments, unlike the expansion of the scalar electrical potential. The first and easiest 
representative of the toroid family is the anapole moment [Jj - a torus with electric currents 
flowing along its meridians. Magnetic multipole moments are defined by transverse currents 
(currents flowing along the torus' parallels for instance). The necessity of the subdivision is 



(*) E-mail: smischenkoOyalioo . com 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



quite obvious since we can decompose the current density into the foUowing terms according 
to a well known theorem of the vector analysis 0: j{r) — grad{r]) + rot{f), where the first 
term corresponds to the toroid moments and the second to the "usual" magnetic multipole 
moments. Spin systems of finite sizes and magnetic molecules (Mni2Ac, V15, Fes, Feio, Feao, 
etc.) in particular are natural objects for application of the multipole expansion technique. 
The magnetic molecules have attracted great attention from the point of view of fundamental 
problems of quantum mechanics in general and the theory of magnetism in particular |3j 
as well as nanotechnology and microelectron applications (as model systems for quantum 
informatics for example) |3l4l5l6i7] in recent years. The molecules are large organic molecular 
complexes with the weight of approximately 10^ atomic mass units. The magnetic ions, such 
as Fe, Mn, V, etc. embedded within the molecule cause their interesting magnetic properties. 
There is a strong exchange interaction between magnetic ions (in the order of 10^ Oe) within 
the molecule. 

Magnetic multipole expansion. - The lowest magnetic multipole moments of spin systems 
(dipole 1, quadrupole 2, octopole 4 and toroid 1) can be defined as follows IHIEIEII: 

Ml = fJ,B^ giSi, {M2)^p ^ l^B^ giinpSia + SifSna), (2) 

i i 

i^3)al3'y = 3^*3 ^ gtStanpri^, (3) 

i 

Mn,m = -f^B^ gi['^i{n"'Cn'^{Oi, fi)) X sl], fi = ^Jlb'^ gi[sifi\, (4) 

i i 

where Mn^m are magnetic multipole moment's components in spherical coordinates {n defines 
the order of a moment: n = 1 - dipole, n = 2 - quadrupole, n = 3 - octopole, eic, m ranges 
from — n to n by 1). It is a well-known fact that each multipole moment interacts with the 
corresponding magnetic field spatial variation. For the sake of brevity we present here the 
formula for the octopole moment only: W = —{l/6){M3)a/3jBjpa, B^^a = V^V^Sq, where 
W is the energy of interaction between an octopole and external magnetic field in Cartesian 
coordinates, and V^g are Laplace operators for the spatial derivatives and Ba stands for 
the three components of the external magnetic field. 

Octopolar electric fields have already been well examined in connection with non-linear op- 
tical properties Quadrupole moment of the rare earth ions, their dependence on external 
magnetic field and connection with magnetic birefringence and magnetostriction have been 
investigated in |12II13II14| . The spin toroid moment Ti has been experimentally found in bulk 
magnets (GaFeOa, Cr203, etc.) recently ^3]. Molecular complexes with the spin density that 
can be described by means of the toroid moment Ti have been analyzed theoretically in '16' 
but an experimental realization of the proposed idea still remains unknown. Quadrupole 
magnetic fields generated by antiferromagnetic crystals have been analyzed in Antifer- 
romagnetic nanoclusters, e.g. Feio,Fe6 can be characterized by the quadrupole moment AI2 
and moments of higher orders. 

Highly symmetrical Keplerate molecule Feao. - However, no one magnetic object has 
been characterized by the octopole moment so far. The first object of this kind is a magnetic 
nanocluster Feao with the chemical formula 

[Mo72Fe30O252(MO2O7(H2O))2(MO2O8H2(H2O))(CH3COO)i2(H2O)9i] *150H2O. 

This molecule has the largest number of magnetic Fe^^ ions among all the known molecular 
magnets synthesized so far The molecule Feao as well as the other nanoclusters discussed 
above occupies an intermediate position among the bulk materials and localized magnetic 
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Fig. 1 



Fig. 2 



Fig. 1 - A schematic picture of the Feao molecule. The three possible spin directions are marked by 
black, gray and white colors. 



Fig. 2 - A model of an elementary octopole consisting of four equal magnetic moments m situated in 
the vertices of the square with a side S. 



ions, that is the reason these objects are called "mesoscopic magnets". Toroid (Ti), dipole 
(All) and quadrupole (M2) moments of this molecule are equal to zero, i.e. the multipole 
expansion starts from the octopole moment. In other words, probably the Fe3o molecule with 
a highly symmetrical spin density is the most symmetrical magnetic body synthesized so far. 

The investigated molecular complex consists of 30 Fe^~^ ions each of which has a spin 
s = 5/2. The magnetic ions are situated in the vertices of the icosidodecahedron, one of 
the regular Archimedean polytopes (see figCJ. The icosidodecahedron has v = 30 vertices, 
f — 32 faces and e = 60 edges in accordance with the so-called Euler theorem [Q, connecting 
a number of vertices, faces and edges of a convex polyhedron of O genus {v + f — e = 2). 
The faces of the icosidodecahedron consist of 20 triangles and 12 pentagons. Neighbouring 
Fe^^ ions interact with each other antiferromagnetically by means of the indirect exchange 
interaction which results in zero total spin of the molecule in the absence of magnetic field . 
The exact ground state energy and the corresponding ground state of the Peso molecule has 
been found in (71 basing on graph-theoretical technique and a special geometrical property of 
the icosidodecahedron, its three-colorability (a graph is called three-colorable if we can color 
all its vertices with three colors and no neighboring vertices with the same color - see fig0. 
In the framework of the proposed model all spins are coplanar and the relative angle between 
nearest-neighbor spins is 120°. In other words, each of the three colors is assigned to one of 
three coplanar spin directions (frgQ. It is important to note that ab-initio calculations of the 
molecule in question have not been performed yet and we have no experimental evidence on 
how the plane coplanar to the spins is aligned - it can be parallel either to one of the pentagons 
or to one of the triangles. However, our calculations show that the result is approximately 
the same in the both cases. The octopole moment components in the spherical system can be 
found from I^J by substituting n = 3; m changes from -3 to 3 by 1, thus yielding seven complex 
components of the moment. Our calculations show that the octopole moment components 
both in Cartesian and spherical systems are finite but we do not present here the particular 
values for the sake of brevity. 

Magnetic field with octopolar symmetry at least is required to measure the octopole mag- 
netic moment components. Such a field can be generated by an elementary magnetic octopole 
presented on fig|21 The field will create a torque on the Keplerate molecule(s) placed in the 
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origin of the coordinate system with the value determined by corresponding octopole moment 
components. The estimation of particular torques is beyond the scope of this article. This 
approach can be used to study a wide number of mesoscopic magnets. Neutron spectroscopy 
can be used to measure the octopole moment as well and further work in this direction is 
required. 

" Deoctopolization" in external magnetic field. - Magnetic properties of the Feao molecule 
and its multipole moments are discussed in the following. However, the dependence of the 
magnetic dipole moment of the magnet in question on external magnetic field has been quan- 
titatively described in 7 in the framework of the classical statistical physics (namely, the 
magnetic moment rises linearly with magnetic field till the saturation at the critical field 
value Be — 17. 7T, what is in perfect agreement with experimental results). The considera- 
tion was based on a study of a three spin system with antiferromagnetic interactions between 
the closest neighbour spins placed in external magnetic field in accordance with the following 
Hamiltonian [18) : 

H = Si S2 + S2 S3 + S3S1 -p{Si+ S2 + S3), (5) 

where (3 = g^iB/2JS is a dimensionless magnetic field, J > stands for antiferromagnetic 
exchange integral and all spins are of the same length S. The Hamiltonian represents a classical 
model, i.e. the projections of the spins are not quantized and can take all intermediate values. 
The exact solution for the ground state energy has been found in jl8| and it was shown that the 
ground state energy depended on the total spin only (or magnetic field, as they are connected 
linearly), namely EasiP) = -3/2 - 0^/2. At the magnetic field value corresponding to /3 = 3 
all the three spins are aligned with the magnetic field. The dimension of the Hilbert space for 
Feao molecule in the framework of the quantum model is 6'^'^ which is far beyond calculation 
possibilities of any modern computer. However, an approximate quantum model based on the 
isotropic nearest-neighbour antiferromagnetic Heisenberg exchange is presented in and 
finally results in appearance of special small quantum steppes on magnetization typical for 
molecular magnets. Analogous consideration would lead to the steppes of the same kind on 
the octopole moment components as well but they would not change the general picture and 
would be an unnecessary complication. This is the reason we use classical model here. 

However, due to the antiferromagnetic interactions the three spin system is highly frus- 
trated and the magnetization process can be carried out in many ways. The three easiest of 
them are: A) the three spins simultaneously leave the plane where they lied before in the case 
of magnetic field directed perpendicularly to the spin plane (figOl a); B) the field is directed 
along one of the spins, then the remaining two spins turn in the field's direction uniformly 
(fig|2| b); C) the field is opposite to one of the spins {Si), then the magnetization process 
involves two phases separated by the Si spin jump at magnetic field (3 — 1 (fig^. 

There is one parameter describing each of the first two magnetization processes: the angle 
of each spin's declination from the initial state i?_l(/3) in case (A); the angle of the ^2 and S3 
spins' declination 'd\\{(3) in case (B) and two parameters in the third case (C) - "daiP) and '0b{(3) 
(figEl- The parameters can be found from minimization of the system's energy in external 
magnetic field. This can be done by projecting the Hamiltonian 101 on subspaces containing 
one or two variables in the first two (A) and (B) and the last case (C), respectively. Due to 
simplicity of the considered cases below we give the answers only: 

i?±(/3) = |/3, i?||,,,(/3) = arcsin^^-^ + ^, (/3) = arcsin ^^-^-^ - ^ , (6) 

where /3 changes from up to its saturation (critical) value 3. Note that the second formula 
describes two parameters - ^?||(/3) and db{P). 
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Fig. 3 - Two considered geometries of the Feso molecule magnetization process: (a) - B is perpendic- 
ular to the spin plane (case A); (b) - _B is parallel to one of the three spins (cases B, C). 



In order to calculate the dependence of high muhipole magnetic moments on the magnetic 
field we have made an assumption that the Feso molecule can be divided into 20 triangles of 
antiferromagnetically interacting spins (O [?]• This assumption can be partially justified by 
the fact that the icosidodecahedron is three-colorable without any visible anisotropy and that 
this model's results are in perfect agreement with experiments. Besides, according to recent 
calculations with the Density Matrix Renormalization Group technique (DMRG) reported 
in j^, the low lying energy levels of the Feao molecule form a rotational band, i.e. they 
depend approximately quadratically on the total spin quantum number S. Actually /j in 
should be multiplied by a factor of 1/2 when writing the Hamiltonian for all 20 triangles 
to take into account that each spin is shared by two neighbouring triangles. Hence, having 
the magnetization parameters © we obtain the spin distribution over the molecule at finite 
magnetic field and can calculate its multipole moments for any field value. As we have already 
noted, the magnetic dipole moment rises linearly with the field. Surprisingly, the quadrupole 
magnetic moment remains equal to zero in all cases (A, B, C) at any fields up to saturation, 
and further. 

The two basic types of the octopole moment components' behavior in the out-of-plane 
case (A) are plotted on figjSl a. As you can see, almost all the components vanish with the 
field. There are only three linearly raising components. It can be understood by a simple 
mathematical consideration. Indeed, at the saturation all the spins are aligned parallel to the 
magnetic field (and Oz axis, see figEl a), i.e their x— and y— components are equal to zero 
as well as 18 octopole magnetic moment components (M3)iij, (M3)2ij, « = l--3,j = 1..3 (jSJ. 




Fig. 4 - A simple scheme of the magnetization process involving the Si spin jump (case C) . 
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Moreover, all components of the kind (M3)3y with i ^ j are equal to zero as the contributions 
from different spins cancel each other out when we perform the summation in 0. The only 
three surviving components are (-^3)311, (M3)322 and (Af3)333 as they are the sums of squared 
coordinates of all spins multiplied by the z component of the saturated spin and it is obvious 
that these three components raise with the field. 

However, the in-plane magnetization cases (B, C) are more interesting. Their general 
geometry is shown on fig|3 b. Note that the magnetic field can have any sign so we cover 
both B and C cases with this picture. Nevertheless, the octopole moment behaves in different 
ways and we will consider these cases separately. Take a look at the fig|Sl b. It differs from 
flgEl ^ with the curvature of M212, M213 components and appearance of dependencies of a new 
type M121, M123. These facts are connected with the way the spin projections change during 
the magnetization process (see © and figEl b). And just like in the (A) case we have three 
linearly growing components - Mm, M122, M133. The third case (C) is the most interesting 
as one of the spins, 5*1 for clarity, abruptly changes its direction for the opposite one at the 
magnetic field /3 = 1 (see figEJ. At the same time, the total spin (and magnetization) do not 
experience a jump or even a kink, what follows from the general solution for the considered 
problem given in jl8|. The octopole moment's components vs. magnetic field are shown on 
figsEla, b. Quite obviously, all the dependencies have jumps at the field /3 = 1 showing the 
abrupt changes of symmetry at this point. 

It is important to point out that the three considered cases do not exhaust all possible ways 
the magnetization can be realized. They rather reveal basic features of the "deoctopolization" 
process and can be a ground for its future classification. 

Conclusion. - The results of the theoretical investigation of the Fe3o molecule let us 
make some important conclusions. First of all it is shown that dipole, toroid and quadrupole 
moments of the Fe3o molecule are equal to zero. Probably the considered system has the 
most symmetric distribution of spin density of all the mesoscopic magnets synthesized so far. 
The magnetic octopole moment is calculated in both Cartesian and spherical coordinates. 
A model of an elementary magnetic octopole as a possible experimental unit for octopole 
magnetic moment component measurement is proposed. Toroid and quadrupole moments are 
shown to remain equal to zero at any magnetic fields while components of the octopole moment 
are changed substantially. Different geometries of the magnetization process are considered 
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and discussed. 
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